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Abstract 

One  candidate  for  the  development  of  a  tunable  laser  in  the  mid-infrared  for  use  as  a 
defensive  countermeasure  against  passive  sensors  exploits  the  nitric  oxide  (NO)  transition  at  5.4 
pm.  However,  the  difficulty  with  this  mechanism  is  that  when  bromine  (Br2)  and  NO  are  present 
in  a  mixture,  nitrosyl  bromide  (BrNO)  is  also  created.  The  absorption  bands  present  in  BrNO 
complicate  the  design  of  such  a  laser.  This  thesis  determines  the  forward  and  reverse  rate 
constants  of  the  reaction  of  Br2  and  NO  to  form  BrNO  following  the  photolysis  of  BrNO  using 
Fourier  Transform  Infrared  absoiption  spectroscopy.  Other  studies  have  been  conducted  to 
determine  both  the  forward  and  reverse  reaction  rate  constants  using  different  methods,  however 
some  of  the  methods  used  create  uncertainty  in  the  results. 

In  this  effort  samples  with  a  total  pressure  of  10.0  torr  were  used  containing  either 
bromine  or  nitric  oxide  in  excess.  The  use  of  Br2  versus  NO  as  the  excess  reagent  has  an  effect 
on  which  mechanism  can  be  used  to  describe  the  reaction,  but  does  not  affect  the  reaction  rate 
constants.  When  NO  is  the  excess  reagent,  the  reaction  obeys  pseudo  first-order  kinetics,  where 
the  reverse  reaction  does  not  significantly  impact  the  overall  reaction.  This  is  not  true  when  Br2 
is  the  excess  reagent.  In  both  cases,  third-order  kinetics  accurately  describe  the  Br2  +  2NO  <-> 
2BrNO  reaction  giving  a  forward  reaction  rate  constant  of  kf  =  1.56  ±  0.20  x  10'38  cm6/molecule2- 
s  at  293  ±  1  K.  The  reverse  rate  constant  was  calculated  as  kr  =  2.29  ±  0.33  x  10'21  cmVmolecule- 
s  and  the  equilibrium  constant  as  Keq  =  171  ±  13  atm"1.  This  result  is  consistent  with  previous 
results. 
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THE  KINETICS  FOLLOWING  PHOTOLYSIS  OF  NITROSYL  BROMIDE 


I.  Introduction 


Problem  Statement 

One  candidate  for  the  development  of  a  tunable  laser  in  the  mid-infrared  for  use  as  a 
defensive  countermeasure  against  passive  sensors  exploits  the  nitric  oxide  (NO)  transition  at  5.4 
pm. 1 1  This  system  utilizes  the  photolysis  of  an  electronically  excited  bromine  atom  (Br*)  which 
collisionally  transfers  energy  to  NOv=2,  the  receptor  molecule.  This  excited  molecule  relaxes 
back  to  the  ground  state  via  radiation,  forming  the  laser.  The  rotational  distribution  within  the 
excited  vibrational  states  provides  a  tunable  range  of  approximately  1  micron.6  The  Br*  in  this 
system  is  produced  by  the  photolysis  of  molecular  bromine  (Br2).n  However,  the  difficulty  with 
this  mechanism  is  that  when  Br2  and  NO  are  present  in  a  mixture,  nitrosyl  bromide  (BrNO)  is  also 
created.  The  absoiption  bands  present  in  BrNO  complicate  the  design  of  such  a  laser.  The 
kinetics  reported  here  help  to  increase  our  understanding  of  the  chemical  reactions  involved  in  the 
production  and  destruction  of  BrNO. 

Additionally,  halogens,  including  Br2,  have  been  linked  to  chemical  reactions  believed  to 
describe  the  ozone  depletion  cycle.  Chlorofluorocarbons  contain  halogens  that  deplete  the  ozone 
through  the  reactions 


X  +  C>3  ->  xo  +  02 

(1) 

XO  +  NO  — >  N02  +  X 

(2) 

o+no2->no+o2 

(3) 

where  X  represents  fluorine,  chlorine,  iodine  or  bromine.  The  net  reaction  reduces  the  ozone 
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concentration. 


O  +  03  ->  20 2  (4) 

The  halogen  atom  is  not  consumed  in  these  reactions  and  may  react  with  nitric  oxide  and  a  third 
body,  Y,  to  form  a  nitrosyl  halide 


X  +  NO  +  Y  — >  XNO  +  Y  (5) 

The  halogen  atoms  also  form  diatomic  molecules  that  react  with  the  NO  again,  forming  additional 
nitrosyl  halides 


X2  +  2NO  <-»  2XNO  (6) 

If  the  halogen  is  present  only  as  X2  or  XNO,  versus  X,  molecules  it  cannot  aid  in  the  depletion  of 
ozone.  A  thorough  understanding  of  the  kinetics  of  the  reactions  is  required  to  determine  the 
effects  of  these  processes.4  This  thesis  addresses  some  of  the  chemical  kinetic  issues  present 
when  a  halogen  reacts  with  nitric  oxide. 

Previous  Results 

Previous  works  studied  the  reaction  of  halogens,  specifically  bromine,  with  nitric  oxide, 
however  some  of  the  methods  used  cause  uncertainty  in  the  results.  Hisatsune  and  Zafonte 
determined  the  forward  and  reverse  reaction  rate  constants  from  the  reaction  of  BrNO,  NO  and 
Br2  by  mixing  NO  and  Br2  and  following  the  time  dependence  of  the  strong  NO  stretch  infrared 
absorption  band  of  the  product  BrNO.8  Using  the  equilibrium  state  of  Br2,  NO  and  BrNO, 
Houel  and  van  den  Bergh  determined  the  UV/visible  absorption  spectrum  of  BrNO.10  Hippier 
studied  the  recombination  of  Br  atoms  after  photolysis  of  less  than  0.3  torr  of  Br2  at  room 
temperature,  in  the  range  of  1-100  atm  of  the  inert  diluent  helium.  He  also  studied  the 
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recombination  of  Br  and  NO  in  the  presence  of  helium.7  Godfrey  examined  the  kinetics  of  BrNO 
formation  and  destruction  using  time-resolved  photolysis  techniques.  The  kinetic  mechanism  and 
rates  for  the  photolysis  of  the  BrNO  mixture  were  also  examined  and  validated  for  total  pressures 
ranging  from  28-111  torr.4  Only  Godfrey  included  the  effect  of  the  loss  of  Br2  to  the  internal 
surfaces  of  the  cell  in  the  calculation  of  the  reaction  rate  constants.  The  previously  determined 
reaction  rates  range  from  1.32  ±  0.14  to  1.68  ±  0.1 1  x  10"18  cm6/molecule2-s  for  kf  and  from  2.09 
±  0.55  to  3.71  x  10'21  cmVmolecule-s  forkr. 
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II.  Theory 


BrNO  Spectrum 

The  fundamental  vibrational  band  centers  of  BrNO  exist  at  1798.7  cm'1  for  the  Vi  band, 
542.0  cm'1  for  the  v2  band,  and  266.4  cm"1  for  the  V3  band.13  The  v,  band  is  primarily  due  to  an 
NO  bond  stretch  while  v2  and  v3  are  predominantly  mixes  of  NBr  stretch  and  angle  bending.  The 
fundamental  band  of  the  two  main  isotopomers  of  BrNO  (79BrNO  and  lS1BrNO)  is  the  strongest 
band  in  the  infrared  spectrum.  The  Vi  band  is  of  parallel  structure  with  only  A-type  transitions. 
The  P  and  R  branches  are  dominant,  but  weaker  Q-branch  features  also  appear  as  a  shoulder  on 
the  P  branch.  The  line  density  of  the  spectrum  is  very  high,  due  to  the  small  values  of  the 
rotational  constants  B  and  C  and  the  presence  of  two  equally  abundant  isotopomers.3 


Wavenumber,  v  (cm1) 


Figure  1.  Composite  2  cm'1  resolution  transmittance  spectrum  of  Vi  vibrational  band  of 
BrNO  measured  at  various  concentrations  of  BrNO. 


4 


Table  1.  Band  centers  (cm'1)  of  BrNO’  13 


79Br14N160 

81Br14N160 

Br15N160 

Br14N180 

Br15N180 

Relative 

Intensity 

Vi 

1798.751 

1798.742 

1768.3 

1751.5 

1719.5 

VVVS 

V2 

542.0 

527.8 

No  Data 

521.5 

VS 

V3 

266.2 

264 

No  Data 

257.0 

S 

Kinetics  Reactions 

BrNO  photodissociates  into  Br2  and  NO  when  it  is  illuminated  with  a  visible  wavelength 
laser  light  source.  After  the  light  is  removed  the  Br2  and  NO  naturally  recombine,  at  a  rate  kf,  to 
form  BrNO,  while  a  backwards  reaction  regenerates  Br2  and  NO,  until  an  equilibrium 
concentration,  Keq,  is  reached,  with  any  excess  Br2  or  NO  remaining  in  the  mixture. 


BrNO  +  hv  — — — >  Br  +  NO 

(7) 

Br  +  BrNO  — — — >  Br2  +  NO 

(8) 

Br ,  +  2 NO  <  >  2 BrNO 

(9) 

The  rate  of  change  of  the  concentration  of  BrNO  over  time  is  described  by 


d[BrNO\ 


=  -kp  [BrNO],  I  -At,  [Br\  [ BrNO \  +  2k f  [M?], 2  [Br2  ],  - 2 kr  [BrNO], 2 


The  calculation  of  the  forward  and  reverse  reaction  rate  constants  in  equation  (9)  can  be 
simplified  by  considering  only  the  data  collected  after  the  light  source  has  been  shut  off,  where 
Ip=0,  and  the  BrNO  mixture  is  returning  to  equilibrium.  The  second  term  in  equation  (10)  is 
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eliminated  because  [Br]  exists  for  less  than  a  microsecond  making  k|  neglible  when  compared  to 
kfand  kr.  This  reaction  is  still  third-order  and  the  rate  of  change  of  BrNO  over  time  is 


d\BrNO\ 

dt 


2t,[Afo],2[a-2],-2t,[awo],! 


<ii) 


Third-order  bond-breaking  reactions  between  stable  molecules  are  quite  rare  and  all 
known  examples  involve  nitric  oxide  as  one  of  the  reactants.14  The  best  understood  reactions  are 
those  involving  chlorine  and  oxygen.  The  collision  of  three  molecules  is  so  unlikely  that  many 
scientists  believe  the  reaction  consists  of  two  steps.  However,  no  experiment  has  ever  provided 
evidence  for  the  existence  of  any  intermediate  complex.8 

The  normalized  concentration  of  BrNO,  which  is  defined  as  x,  can  be  calculated  from  the 
FTIR  spectroscopy  data  collected  for  each  sample  by  measuring  at  time,  t,  and  at  equilibrium,  the 
area  under  the  curve  of  the  BrNO  (v=l)  band. 


[BrNO], 
[ BrNO]eq 


(12) 


The  Beer-Lambert  Law  describes  the  absorption,  A,  of  radiation  by  molecules  after 
passing  through  a  length,  l,  of  gas,  such  as  during  FTIR  spectroscopy,  as 


A  =  log 


10 


V7o  J 


=  -eel 


or 


A  =  log 


f 


10 


V70  J 


-  a  BrNO  J 
ln(10) 


(13) 


6 


where  the  I  is  the  intensity  of  the  beam  passing  through  the  filled  cell,  Io  is  the  intensity  detected 
when  the  cell  is  empty,  c  is  the  concentration  in  mol/liter,  [BrNO]  is  the  concentration  of  BrNO  in 
molecules/cm\  s  is  the  absoiption  coefficient  in  liter/mole-cm,  a  is  the  absorption  cross  section 
in  cm2,  and  /  is  in  centimeters.9  Since  a  and  /  are  constant,  the  Beer-Lambert  Law  can  be  used  to 
determine  that  the  ratio  of  areas  under  the  absoiption  or  transmittance  curve  for  a  specific  peak 
taken  at  two  different  times  during  the  reaction  is  equal  to  the  ratio  of  the  concentrations  of  the 
sample  at  those  times. 


Figure  2.  Representative  sample  of  measured  2  cm'1  resolution  FTIR  transmittance  data 
at  time  =  0,  time  =  7.5  min,  and  t  =  63  min  (equilibrium).  P0(Br2)  =  0.25  torr, 
P0(NO)  =  9.75  torr,  P0(BrNO)  =  0  torr. 
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The  molecular  concentrations  for  Br2,  NO,  and  BrNO  at  equilibrium  are 


(14) 

[NOl^[NO\-2n,t 

(15) 

[BrNO\,  =[BrNO\+2n„ 

(16) 

where  neq  is  a  constant  used  to  determine  the  equilibrium  concentration  of  each  species  in  the 
sample.  By  substituting  equations  (14),  (15),  and  (16)  into  the  equation  for  the  equilibrium 
constant, 


k,  [BrNOl, 

K  [mT>T, 


a  value  for  neq  can  be  found  from  the  initial  concentrations  of  Br2,  NO,  and  BrNO.  Since 
[BrNO]o  =  0,  the  value  for  neq  is  calculated  from 


Additionally,  the  concentrations  of  Br2,  NO,  and  BrNO  at  time,  t,  in  terms  of  x  are 


[BrNO],  =x[BrNO\t 

(19) 

[NO\=[NOl,+[BrNOlq(\-x) 

(20) 

[B<-A=[BrA,,+l\BrNO\q(\-x) 

(21) 
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Substituting  these  values  into  equation  (11)  gives 


dx 

dt 


2  kf[BrNO]-e 


[NO],, 


[■ BrNO\,q 


+  (!-*) 


-  1  ll-.v)|-2/,.  [BrNO], 
[BrNO\,  2'  ') 


X 


(22) 


By  letting 


[NO\, 

[BrNO\, 

(23) 

K,=2  k,[BrNOl, 

(24) 

K2=2  t,[MO], 

(25) 

the  differential  equation  becomes 


dx 

dt 


K1(a  +  (l-Jc))J 


v 


K,  1  (,  0 

Kj  a  2  j 


-K2x 


(26) 


The  solution  to  the  third-order  problem,  which  was  found  using  Mathematica,  is 


/4 


(27) 


where 


f  =  V-  K2(l  +  «)2(k2(ct  -  l)2  -  2K,ar3 ) 
f2  =  -2(k2  -K2cr  +  K^3 )ArcTan(K,a2{x^Kl^2~l)>) 
/3=/1(-ln((x-l)2)+ln(K1a2(l  -x  +  af  +2K2(l  +  a)(l  +  x(a-l)  +  a))) 
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/4=/,(Kiar3+4K2(l  +  a)) 


Figure  3  shows  a  plot  of  the  third-order  solution.  The  values  of  a,  Ki,  and  K2  determine  the 
shape  of  the  curve.  All  three  values  must  be  positive,  as  defined  by  equations  (23),  (24),  and 
(25),  but  regardless  of  a,  Ki,  and  K2,  the  solution  has  an  asymptote  at  x=l. 


Normalized  BrNO  concentration,  x  (unitless) 

Figure  3.  Solution  to  third-order  kinetics  problem.  Plotted  for  the  case  of 

P0(Br2)  =  0.3  torr,  P0(NO)  =  9.7  torr,  a  =  16.15,  K[  =  1.13E-05  sec'1,  K2  = 
8.56E-05  sec'1. 

Values  for  kf,  kr,  and  Keq  can  be  calculated  by  using  the  estimated  values  of  a,  Ki  and  K2. 
Section  IV  discusses  how  the  solution  was  fit  to  the  collected  spectrometric  reaction  rate  data. 
However,  when  this  solution  is  fit  to  the  collected  data,  the  computed  values  for  a,  Ki,  and  K2 
have  error  values  over  100%  larger  than  the  actual  value  of  the  parameter.  This  error  results 
because  Tablecurve,  the  curve  fitting  program  used  in  this  effort,  has  too  many  fit  parameters  to 
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be  certain  of  the  value  of  any  of  them.  The  number  of  fit  parameters  in  equation  (27)  can  be 
reduced,  consequently  reducing  the  error,  by  letting 


K,=K2*Kj*-JVO],t  (28) 

and  computing  a  from  [NO]0  and  [Br2]o 

g.  fr°l»  ,29, 

-  [BrNO\,  ~  P9) 

where  neq  is  a  constant  determined  from  Keq,  [NO]oand  [Br2]o  by  equation  (18). 

When  an  excess  of  NO  is  present,  the  reaction  of  Br2  +  2NO  •*->  2BrNO  can  be 
approximated  as  a  first-order  reaction,  a  gets  large,  Kia2  becomes  much  larger  than  K2,  and 
equation  (26)  reduces  to 


dx 

dt 


Kja2 


ri_ 

u 


(30) 


The  solution  for  this  equation  is 


x  =  l  +  Cte~*At  (31) 

where  A  =  K^or.  Values  for  kf  and  kr  can  be  calculated  by  using  the  fitted  value  of  A.  Section  IV 
discusses  how  the  solution  was  fit  to  the  data  collected  in  which  there  was  an  excess  of  NO.  Q  is 
calculated  from  the  initial  conditions  present  in  the  data. 

Appendix  A  provides  the  complete  calculation  for  determining  the  kinetics  of  the 
Br2  +  2NO  <-*■  2BrNO  reaction. 
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III.  Experimental  Procedure 


Experimental  Set  Up 

The  experiment  was  conducted  at  10.0  torr  total  pressure  in  a  10  cm  glass  cell  with 
calcium  fluoride  (CaF2)  windows.  A  gas  handling  vacuum  system  was  used  to  fill  the  glass  cell 
with  a  mixture  of  Br2  (99.5%  pure,  Spectrum  Chemical  Manufacturing)  and  NO  (99%  pure, 
Matheson  Gas)  at  room  temperature  (T=293  K)  with  samples  containing  either  an  excess  of  Br2 
or  NO.  The  evacuated  glass  cell  was  first  backfilled  with  helium  to  eliminate  any  residual 
chemicals  that  might  be  in  the  cell  from  previous  samples.  Additionally,  for  several  samples,  the 
cell  was  cleaned  with  methanol  to  eliminate  the  buildup  of  Br2  on  the  walls  of  the  cell.  A  vacuum 
was  pulled  on  the  system  until  the  total  pressure  was  10  millitorr.  After  the  entire  system  was 
evacuated,  the  glass  cell  was  filled  with  the  substance  that  had  the  lower  partial  pressure.  The 
glass  cell  was  isolated,  the  vacuum  system  was  again  evacuated,  and  the  substance  with  the 
higher  partial  pressure  was  added  to  the  cell.  Each  sample  had  a  total  pressure  of  10.0  torr,  but 
the  partial  pressures  of  Br2  and  NO  were  varied  to  demonstrate  that  the  initial  composition  of  the 
mixture  did  not  have  an  effect  on  the  rate  constants. 

The  first  four  samples  were  filled  at  room  temperature  with  1.0  torr  of  NO  and  9.0  torr  of 
Br2  and  the  remaining  samples  were  filled  with  9.0  torr  of  NO  and  1.0  torr  of  Br2,  but  because  of 
the  tendency  of  bromine  to  adhere  to  the  internal  surfaces  of  the  glass  cell  and  CaF2  windows,  the 
actual  ratio  of  bromine  and  NO  in  the  cell  at  the  beginning  of  data  collection  is  not  the  same  as 
the  original  ratio  used  when  mixing  the  samples.  The  amount  of  bromine  lost  is  difficult  to 
quantify  because  the  loss  rate  depends  on  the  quantity  of  bromine  in  the  cell  and  is  different  each 
time  the  cell  is  filled.  The  bromine  begins  to  stick  to  the  walls  almost  immediately  making  it 
nearly  impossible  to  prevent  this  phenomenon.  Instead  of  trying  to  prevent  the  loss  of  Br2  to  the 
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walls  of  the  cell,  the  ratio  of  Br2  to  NO  at  the  start  of  data  collection  for  each  sample  was 
determined  from  the  data  collected. 


experimental 


Figure  4.  Gas  handing  vacuum  system  used  to  fill  the  glass  cell. 

Only  one  of  the  previous  experiments  for  determining  the  forward  and  reverse  rate 
constants  of  the  Br2  +  2NO  <->  2BrNO  reaction  addressed  the  loss  of  bromine  in  the  sample  to  the 
walls.  Flowever,  no  viable  solution  to  the  problem  was  found.  The  loss  of  bromine  to  the  walls 
introduces  more  error  into  the  determination  of  the  rates,  but  in  this  study  the  error  is  accounted 
for  by  the  extrapolation  of  the  pressure  ratio  of  Br2  to  NO  at  the  start  of  data  collection.  Section 
IV  explains  in  more  detail  how  the  pressure  ratio  at  the  start  of  data  collection  is  found  for  each 
sample. 
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To  demonstrate  the  tendency  of  bromine  to  react  with  the  internal  surfaces  of  the  cell, 
pressure  data  was  recorded  with  a  100-torr  Baratron  pressure  manometer  from  two  separate 
samples  containing  1.0  torr  of  bromine  in  a  glass  cell  over  a  period  of  30  minutes.  The  first 
sample  measured  the  pressure  drop  in  both  the  glass  cell  and  the  quarter-inch  stainless  steel 
tubing  of  the  gas  handling  system  while  the  second  sample  measured  the  pressure  drop  in  only  the 
tubing.  The  pressure  drop  in  only  the  glass  cell  was  found  by  taking  the  difference  between  the 
two  samples  because  the  setup  of  the  gas  handling  system  did  not  allow  for  the  pressure  drop  in 
only  the  glass  cell  to  be  measured  directly.  This  test  was  performed  simply  to  show  that  wall 
effects  exist  within  the  sample  cells  used  in  this  experiment.  More  data  is  needed  to  determine 
the  rate  at  which  the  bromine  pressure  decreases  in  the  cell.  The  loss  of  bromine  to  the  internal 
surfaces  of  the  cell  explains  why  the  initial  bromine  pressure  in  the  cell  during  data  collection  is 
lower  than  the  pressure  of  bromine  originally  added  to  the  cell. 


0.80  -| - T - T - 1 - 1 

0  500  1000  1500  2000 

Time,  t  (sec) 

Figure  5.  Pressure  decrease  over  time  due  to  loss  of  Br2  to  the  internal  surfaces  of  the 
system.  Samples  contained  1.0  torr  of  Br2  initially. 
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Photodissociation  and  Reformation  of  BrNO 

Immediately  after  the  cell  was  filled,  the  mixture  of  Br2  and  NO  began  reacting  to  form 
BrNO  as  described  in  equation  (9).  The  Spectra  Physics  Millenia  Vis  J,  a  CW  diode  pumped 
frequency  doubled  Nd:YV04  laser  operating  at  532  nm,  was  used  to  illuminate  the  cell,  causing 
the  photodissociation  of  BrNO  to  form  Br2  and  NO  as  shown  in  equations  (7)  and  (8).  The  laser 
output  was  set  to  4  W  to  maximize  beam  quality.  A  neutral  density  filter  wheel  was  used  in 
conjunction  with  two  cylindrical  lenses,  which  broadened  the  beam,  to  decrease  the  irradiance  on 
the  cell  to  0.1  ±  0.05  mW/cm2  for  samples  with  an  excess  of  NO  and  25  ±  10  mW/cni2  for 
samples  with  an  excess  of  Br2.  The  expanded  beam  was  necessary  to  ensure  all  the  BrNO  in  the 
cell  was  photodissociated.  After  the  mixture  had  reacted  with  the  light  to  form  only  Br2  and  NO, 
the  light  was  removed  from  the  sample  and  the  mixture  was  allowed  to  return  to  equilibrium. 


Figure  6.  Apparatus  used  to  collect  kinetics  data. 
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Collection  of  Spectra 

The  spectra  used  in  this  experiment  were  obtained  using  the  Bomem  MB  157  FT-IR 
Spectrometer,  which  utilizes  a  modulated  broad  band  IR  source  and  DTGS  detector.  The  infrared 
beam  is  modulated  before  passing  through  the  sample.  An  interferometer  causes  each  IR 
frequency  to  be  modulated  with  a  unique  frequency  of  modulation.  After  the  IR  beam  has  passed 
through  the  sample,  the  intensity  is  detected  and  the  frequencies  are  demodulated  via  a  Fourier 
Transform.  The  absorbance  or  transmittance  spectra  are  obtained  after  calculations  are  executed 
to  combine  the  pre-recorded  reference  spectrum  and  the  sample  spectrum. 


Figure  7.  Simplified  schematic  of  the  MB  157  FTIR  Spectrometer. 

The  heart  of  the  FTIR  spectrometer  is  the  Michelson  interferometer.  The  light  from  the 
IR  source  enters  the  interferometer  and  encounters  a  beamsplitter  that  reflects  one  beam  towards 
one  cube  comer  mirror  set,  which  reflects  it  back  towards  the  beamsplitter,  and  transmits  the 
other  beam  towards  the  other  cube  comer  mirror  set,  which  also  reflects  it  back  to  the 
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beamsplitter.  As  both  mirrors  move,  the  optical  path  length  difference  between  the  two  beams 
continuously  changes.  As  a  moving  mirror  is  scanned,  the  two  returned  beams  interfere  with 
different  phases,  which  creates  intensity  variations.  At  a  given  optical  path  length  difference,  the 
interference  is  constructive  at  some  frequencies  and  destructive  at  others.  Because  the  optical 
path  length  difference  is  constantly  changing,  the  various  frequencies  present  in  the  beam  are 
modulated  are  different  rates. 

After  leaving  the  interferometer,  the  modulated  light  passes  through  the  sample  where 
some  of  the  IR  radiation  is  absorbed  at  specific  frequencies.  The  remaining  IR  radiation  reaches 
the  DTGS  detector,  which  converts  it  to  an  electrical  signal  that  is  then  recorded  as  an 
interferogram.  The  spectrum  is  obtained  from  the  interferogram  using  a  Fourier  Transform.2 

Safety  Precautions 

The  primary  hazards  associated  with  this  experiment  were  the  laser  and  the  reactants  used 
(Br2  and  NO).  A  Class  IV  laser  was  operated  at  4  W  to  illuminate  the  sample  of  BrNO  during 
this  experiment.  Permanent  eye  damage  can  result  from  unprotected  exposure  to  laser  radiation. 
Appropriate  laser  safety  glasses  were  worn  and  all  jewelry  was  removed  during  the  use  of  the 
laser.  Additionally,  bromine  and  nitric  oxide  are  serious  health  hazards.  Bromine,  a  dark  red 
liquid  at  room  temperature,  volatizes  rapidly  to  form  a  red  vapor  with  a  strong,  disagreeable  odor. 
The  vapor  is  very  irritating  to  the  eyes  and  throat;  the  liquid  causes  painful  sores  when  spilled  on 
the  skin.  Nitric  oxide  is  a  colorless,  poisonous  gas  that  reacts  rapidly  with  oxygen  to  form  N02,  a 
mildly  poisonous  brown  vapor  that  is  one  of  the  principle  components  of  smog.17  NO  gas 
detectors  monitored  for  leaks  around  the  area  containing  the  gas  handling  system.  Great  care  was 
taken  during  this  experiment  to  avoid  direct  exposure  to  any  of  these  chemicals. 
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IV.  Results 


Third-Order  Solution 

A  2  cm'1  resolution  FTIR  transmittance  spectrum  was  collected  approximately  every  six 
seconds  for  every  sample  cell  from  the  moment  the  laser  was  turned  on  and  began  reacting  with 
the  sample  until  the  sample  reached  equilibrium  after  the  laser  was  removed.  The  spectra 
collected  while  the  laser  was  illuminating  the  sample  were  used  to  verify  that  all  the  BrNO  had 
been  dissociated  into  Br2  and  NO.  The  spectra  collected  after  the  laser  was  turned  off  were  used 

to  determine  the  forward  and  reverse  reaction  rate  constants  for  Br2  +  2NO  < — — — >  2BrNO  . 


Figure  8.  Example  of  area  measured  “under”  the  BrNO  v=l  line  by  Bomem 
GRAMS/32.  Peq(BrNO)  =  0.565  ton- 
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At  equilibrium,  a  very  strong  v=l  line  of  BrNO  is  observed  in  the  FTIR  spectra  collected. 
Using  Bomem  GRAMS/32,  the  area  under  the  BrNO  line  could  be  calculated  for  all  the  spectra 
acquired  from  each  sample.  These  areas  were  normalized  with  the  area  under  the  peak  when  the 
sample  was  at  equilibrium.  As  explained  in  Section  II  the  normalized  value  of  the  area  under  the 
BrNO  peak  is  equal  to  the  normalized  concentration  of  BrNO,  x,  which  is  needed  to  fit  the  data  to 
the  third-order  solution  of  the  reaction.  However,  the  solution  in  equation  (27)  for  the  third-order 
kinetics  problem  cannot  be  solved  for  x  in  terms  of  t,  but  instead  the  data  points  must  be  reversed 
to  allow  equation  (27)  to  be  fit  to  the  data.  The  error  values  of  the  normalized  BrNO 
concentration  for  each  of  the  samples  were  estimated  using  the  standard  deviation  of  the  area 
measurements. 
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Figure  9.  Example  of  normalized  BrNO  concentration,  x,  as  a  function  of  time,  t.  The 
laser  is  turned  off  at  time=0  sec  and  the  sample  is  allowed  to  reach 
equilibrium.  Po(Br2)  =  0.3  torr,  P0(NO)  =  9.7  ton- 
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4000 


Figure  10.  Example  of  third-order  solution  fit  (r2  =  0.952)  to  observed  data  with  K2  as 
the  only  fit  parameter.  P0(Br2)  =  0.28  torr,  P0(NO)  =  9.72  torr 


Using  Tablecurve,  a  curve  fitting  program,  a  value  for  K2  was  found  for  each  of  the 
samples.  The  value  for  K2  is  dependent  on  the  value  of  a,  which  in  turn  is  dependent  on  [Br2]0 
and  [NO]o  according  to  equation  (29).  The  values  for  [Br2]0  and  [NO]o  needed  to  be  established 
from  the  data  collected  because  of  the  loss  of  Br2  to  the  walls  of  the  cell.  Using  a  trial  and  error 
method,  the  total  pressure  in  the  cell  was  held  constant  at  10.0  torr,  while  the  initial  fraction  of 
NO  (f(NO)o)  and  the  initial  fraction  of  Br2  (f(Br2)0)  were  varied.  The  values  used  for  [Br2]0  and 
[NO]0  were  the  values  which  gave  the  best  curve  fit  to  the  data  when  calculating  K2. 
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Table  2.  Example  of  calculation  of  [Br2]o  and  [NO]0  from  the  curve  fit  data. 
PTo  =  10.0  torr,  T=293  K,  K^q  =  6.71E-18  cm3/molecule 


f(NO)0 

(unitless) 

f(Br2)0 

(unitless) 

[BrNO\q 

(unitless) 

Keq[BrNO]eq 

(unitless) 

K, 

(sec-1) 

k2 

(sec-1) 

K2  Error 

(sec-1) 

r2 

0.900 

0.100 

4.67 

0.351 

— 

won't  fit 

— 

— 

0.910 

0.090 

5.21 

0.324 

— 

won't  fit 

— 

— 

0.920 

0.080 

5.88 

0.296 

— 

won't  fit 

— 

— 

0.930 

0.070 

6.74 

0.266 

— 

won't  fit 

— 

— 

0.940 

0.060 

7.91 

0.233 

— 

won't  fit 

— 

— 

0.950 

0.050 

9.55 

0.199 

2.18E-05 

1.10E-04 

5.33E-06 

0.744 

0.960 

0.040 

12.02 

0.163 

1.49E-05 

9.13E-05 

2.82E-06 

0.886 

0.965 

0.035 

13.79 

0.144 

1.20E-05 

8.32E-05 

2.00E-06 

0.928 

0.970 

0.030 

16.16 

0.125 

9.44E-06 

7.55E-05 

1.50E-06 

0.949 

0.972 

0.028 

17.34 

0.117 

8.48E-06 

7.24E-05 

1.39E-06 

0.952 

0.975 

0.025 

19.48 

0.105 

7.12E-06 

6.77E-05 

1.32E-06 

0.950 

0.977 

0.023 

21.22 

0.097 

6.26E-06 

6.44E-05 

1.32E-06 

0.945 

0.980 

0.020 

24.46 

0.085 

5.05E-06 

5.93E-05 

1.36E-06 

0.930 

0.990 

0.010 

49.43 

0.043 

1.67E-06 

3.84E-05 

1.38E-06 

0.826 

Table  3.  Summary  of  partial  pressures  at  the  start  of  data  collection.  Values 
extrapolated  from  the  measured  FTIR  data. 


Sample 

Total  Pressure 
(torr) 

P(NO)0 

(torr) 

P(Br2)0 

(torr) 

1 

10.0 

0.93 

9.07 

2 

10.0 

0.92 

9.08 

3 

10.0 

0.92 

9.08 

4 

10.0 

0.92 

9.08 

5 

10.0 

9.75 

0.25 

6 

10.0 

9.70 

0.30 

7 

10.0 

9.75 

0.25 

8 

10.0 

9.77 

0.23 

9 

10.0 

9.75 

0.25 

10 

10.0 

9.70 

0.30 

11 

10.0 

9.70 

0.30 

12 

10.0 

9.75 

0.25 

13 

10.0 

9.65 

0.35 

14 

10.0 

9.65 

0.35 

15 

10.0 

9.65 

0.35 
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After  determining  the  pressure  values  of  Br2  and  NO  at  the  start  of  data  collection  for 
each  sample,  the  values  of  kf  and  kr  were  calculated  using  the  initial  concentrations  of  Br2  and 
NO,  which  are  easily  calculated  from  the  initial  pressures  using  the  ideal  gas  law,  the  equilibrium 
constant,  and  the  value  of  K2  found  by  fitting  the  third-order  solution  to  the  data.  The  mean  value 


±  la,  calculated  from  the  data  in  Table  4,  of  kf  is  1.55  ±  0.19  x  10"38  cm6/molecule2-s  and  of  kr  is 


2.31  ±  0.28  x  10'21  cmVmolecule-s. 

Table  4.  Calculated  kf  and  kr  values  using  data  fit  to  third-order  solution  with  K2  as  only 
fit  parameter. 


Po(Br2) 

(torr) 

Po(NO) 

(torr) 

kf 

(10‘38  cm6/molecule2-s) 

kr 

(1021  cm3/molecule-s) 

9.07 

0.93 

1.395 

2.079 

9.08 

0.92 

1.398 

2.084 

9.08 

0.92 

1.342 

2.000 

9.08 

0.92 

1.236 

1.843 

0.25 

9.75 

1.411 

2.104 

0.30 

9.70 

1.472 

2.194 

0.25 

9.75 

1.605 

2.393 

0.23 

9.77 

1.617 

2.411 

0.28 

9.72 

1.391 

2.074 

0.30 

9.70 

1.595 

2.377 

0.30 

9.70 

1.604 

2.391 

0.25 

9.75 

1.735 

2.586 

0.35 

9.65 

1.693 

2.523 

0.35 

9.65 

1.885 

2.810 

0.35 

9.65 

1.836 

2.736 

The  values  calculated  for  kf  and  kr  using  the  third-order  solution  with  only  K2  as  a  fit 
parameter  are  refined  by  fitting  both  Ki  and  K2  to  the  third-order  solution  with  a  fixed  value  for 
a.  Keeping  a  fixed  holds  the  error  in  the  computed  values  of  Ki  and  K2  to  less  than  10%  of  the 
value.  However,  any  attempt  to  fit  all  three  variables,  a,  Kls  and  K2,  to  the  collected  kinetics 
data,  gives  error  values  of  greater  than  100%  for  all  of  the  fit  parameters,  causing  the  fit 
parameters  to  be  essentially  useless.  Therefore,  the  fit  of  K!  and  K2  to  the  kinetics  data  gives  the 
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most  reliable  values  for  kf  and  kr.  Additionally,  the  value  of  Keq  can  be  calculated  from  the  values 
determined  for  kf  and  kr.  The  values  of  the  rate  constants  and  the  equilibrium  constant 
determined  in  this  experiment  are  listed  in  Table  6  along  with  the  results  from  previous  efforts. 
The  data  used  to  determine  kf,  kr,  and  Keq  is  shown  in  Table  5. 


Table  5.  Calculated  Keq,  kf  and  kr  values  using  data  fit  to  third-order  solution  with  K, 
and  K2  as  fit  parameters. 


P0(Br2) 

(torr) 

Po(NO) 

(torr) 

kf 

(10'38  cm6/molecule2-s) 

kr 

(10‘21  cmVmolecule-s) 

Keq 

(atm'1) 

9.07 

0.93 

1.483 

2.217 

168 

9.08 

0.92 

1.372 

2.044 

168 

9.08 

0.92 

1.326 

1.976 

168 

9.08 

0.92 

1.184 

1.763 

168 

0.25 

9.75 

1.394 

2.152 

162 

0.30 

9.70 

1.542 

2.011 

192 

0.25 

9.75 

1.614 

2.361 

171 

0.23 

9.77 

1.628 

2.372 

172 

0.28 

9.72 

1.409 

2.023 

174 

0.30 

9.70 

1.628 

2.298 

177 

0.30 

9.70 

1.724 

2.151 

201 

0.25 

9.75 

1.647 

2.826 

146 

0.35 

9.65 

1.664 

2.618 

159 

0.35 

9.65 

1.867 

2.865 

163 

0.35 

9.65 

1.861 

2.642 

176 

Table  6.  Values  of  kf,  kr,  and  Keq  obtained  in  this  and  previous  efforts.  The  quantities 
listed  for  this  study  are  the  means  ±  la  determined  from  the  data  in  Table  5. 


kf(10'38 

cm6/molecule2-s) 

kr(10'21 

cm3/molecule-s) 

Keq  (atm'1) 

T(K) 

This  Effort 

1.56  ±0.20 

2.29  ±0.33 

171  ±  13 

293  ±  1 

Godfrey4 

1.40  ±0.18 

2.09  ±0.55 

168  ±23 

293  ±  1 

Hippier7 

1.68  ±  0.1 1 

— 

— 

298 

Houel10 

1.60  ±0.16 

2.67 

150  ±48 

303  ±  0.5 

Hisatsune8 

1.32  ±0.14 

3.71 

89 

303 
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First-Order  Solution 


The  forward  and  reverse  reaction  rate  constants  could  also  be  calculated  using  the  first- 
order  solution,  equation  (31),  for  samples  that  contained  an  excess  of  NO.  The  parameters  A, 
which  equals  K,a2,  and  Ci  were  fit  to  the  data.  Mathematically,  at  the  point  x=0,  t=0,  Ci  is  equal 
to  -1.  In  this  experiment,  the  data  determined  that  Q  was  approximately  -1,  showing  that  the  data 
was  very  close  to  what  was  expected.  Due  to  system  noise,  the  concentration  of  BrNO  is  never 
actually  measured  as  zero,  even  after  all  the  BrNO  has  reacted  to  form  NO  and  Br2. 

The  variable  a  was  calculated  using  the  values  of  [Br2]0  and  [NO]o  determined  while 
fitting  the  third-order  solution  to  the  data.  kf  was  then  calculated  from  the  fit  parameter  A  and  kr 
was  found  using  kf  and  the  equilibrium  constant.  The  mean  value  ±  1  a,  calculated  from  the  data 
in  Table  7  of  kf  is  1.91  ±  0.69  x  10'38  cm6/molecule2-s  and  of  kr  is  2.84  ±  1.03  x  10‘21 
cmVmolecule-s. 


Table  7.  Calculated  kf  and  kr  values  using  data  fit  to  first-order  solution  with  AMCjOi2 
and  Ci  as  fit  parameters. 


P0(Br2)  (torr) 

Po(NO)  (torr) 

kf(10‘38  cm6/molecule2-s) 

kr(10"21  cmVmolecule-s) 

0.25 

9.75 

1.42E-38 

2.12E-21 

0.30 

9.70 

1.62E-38 

2.42E-21 

0.25 

9.75 

1.42E-38 

2.11E-21 

0.23 

9.77 

1.30E-38 

1.94E-21 

0.25 

9.75 

1.43E-38 

2.13E-21 

0.30 

9.70 

1.67E-38 

2.49E-21 

0.30 

9.70 

1.63E-38 

2.43E-21 

0.25 

9.75 

1.60E-38 

2.39E-21 

0.35 

9.65 

3.00E-38 

4.48E-21 

0.35 

9.65 

2.92E-38 

4.36E-21 

0.35 

9.65 

2.96E-38 

4.42E-21 
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Figure  11.  Example  of  first-order  solution  fit  (r2  =  0.997)  to  observed  data  with 

A=K!a2  and  Ci  as  fit  parameters.  Po(Br2)  =  0.25  torr,  Po(NO)  =  9.75  torr 


Solution  Near  Zero  BrNO  Concentration 

The  original  differential  equation,  equation  (26),  becomes  equation  (32)  as  x  approaches 

zero. 


dx 


dt 


x=0 


==Kj(a  +  l)- 


K. 


vK,cr 


1 

■  +  — 


(32) 


Table  8  illustrates  the  difference  between  the  theoretical  slope  calculated  with  this 
equation  and  the  slope  measured  from  the  data.  The  slope  of  each  sample  was  measured  near  x=0 
by  computing  the  linear  regression  of  the  first  50  data  samples,  which  occurs  when  time  is  less 
than  286  seconds. 
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Figure  12.  Example  of  concentration  versus  time  near  [BrNO]t  equals  zero. 

Po(Br2)  =  0.30  torr,  P0(NO)  =  9.70  torr.  Measured  slope  near  concentration 
equals  zero  is  0.0013  sec'1.  Calculated  slope  equals  0.0015  sec'1. 


Table  8.  Measured  versus  calculated  slope  near  [BrNO]t  =  0 


Sample 

a 

(unitless) 

Im  (sec'1) 

K2  (sec'1) 

j  measured 

ax 

dt  ,-=o 

(sec"1) 

7  theoretical 

ax 

^  x=0 
(sec'1) 

1 

0.72 

9.37E-06 

7.87E-05 

0.00039 

0.00047 

2 

0.72 

8.55E-06 

7.22E-05 

0.00028 

0.00043 

3 

0.72 

8.30E-06 

6.99E-05 

0.00033 

0.00041 

4 

0.72 

7.37E-06 

6.23E-05 

0.00032 

0.00037 

5 

19.48 

7.11E-06 

6.75E-05 

0.00112 

0.00156 

6 

16.16 

9.36E-06 

7.49E-05 

0.00128 

0.00146 

7 

19.48 

7.80E-06 

7.41E-05 

0.00111 

0.00172 

8 

21.21 

6.68E-06 

6.88E-05 

0.00093 

0.00172 

9 

19.48 

7.85E-06 

7.46E-05 

0.00117 

0.00173 

10 

16.16 

1.07E-05 

8.56E-05 

0.00125 

0.00167 

11 

16.16 

1.00E-05 

8.01E-05 

0.00133 

0.00156 

12 

19.48 

7.85E-06 

7.46E-05 

0.00121 

0.00173 

13 

13.79 

1.62E-05 

1.13E-04 

0.00228 

0.00190 

14 

13.79 

1.78E-05 

1.23E-04 

0.00218 

0.00208 

15 

13.79 

1.64E-05 

1.14E-04 

0.00223 

0.00192 
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Ideally  the  theoretical  slope  would  equal  the  measured  slope,  but  due  to  errors  in 
calculating  both  the  measured  and  theoretical  slopes,  the  values  are  slightly  different.  Error  is 
introduced  into  the  calculation  of  the  measured  slope  from  the  choice  of  endpoints  and  the 
variability  of  the  measured  data.  Due  to  noise  in  system,  the  normalized  concentration  is  never 
measured  as  zero,  which  causes  most  of  the  error  in  the  measured  slope.  The  calculated 
theoretical  slope  contains  the  statistical  error  in  K,  and  K2. 

Error 

Error  in  the  rates  is  mostly  due  to  the  uncertainty  in  the  initial  pressure  measurements. 
The  error  in  the  pressure  of  Br2  and  NO  added  to  the  cell  is  ±  0.01  torr,  due  to  the  limitations  of 
the  manometer  used  to  monitor  the  pressure.  The  error  in  the  pressures  of  Br2  and  NO  at  the  start 
of  data  collection  is  larger,  however,  because  the  values  had  to  be  extrapolated  from  the  collected 
kinetics  data  due  to  the  loss  of  Br2  to  the  walls  of  the  cell.  The  error  in  determining  the  initial 
values  of  Br2  and  NO  is  ±  0.1  torr. 

Temperature  variations  also  change  the  pressure  in  the  cell.  The  temperature  inside  the 
cell  was  not  measured,  but  the  relatively  low  heat  capacity  of  the  gas  mixture  permits  the 
assumption  of  thermal  equilibrium  inside  the  cell.  The  heat  capacity  of  the  glass  cell  also  allows 
the  assumption  that  sudden  temperature  variations  within  the  cell  did  not  occur.  Regardless,  a  ±  1 
K  temperature  variation  at  293  K  changes  the  pressure  of  a  10  torr  ideal  gas  by  ±  0.034  torr, 
which  is  small  when  compared  to  the  ±0.1  torr  error  in  the  extrapolated  initial  pressure  values. 

The  statistical  error  in  fitting  the  parameters,  K2  and  K1?  to  the  measured  data  was  another 
source  of  error.  The  amount  of  error  in  the  fit  parameters  to  the  third-order  solution  was  small 
compared  to  the  actual  value  of  the  parameter,  at  less  than  10%  of  the  value,  whereas  the  error  in 
the  fit  parameters  to  the  first-order  solution  was  on  average  over  20%  of  the  value.  Additionally, 
Tablecurve  calculated  how  well  the  theoretical  solution,  containing  the  computed  values  of  K2 
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and  Ki,  fit  to  the  measured  data  using  both  the  Linear  Least  Squares  Method  and  the  Levenberg- 
Marquardt  (nonlinear  least  squares)  Method.  For  every  sample  an  r2  value  was  calculated  using 
both  methods,  where  r2  =  1  indicates  that  the  equation  fits  the  data  perfectly.  In  this  experiment 
both  methods  gave  identical  r2  values  to  two  significant  figures.  For  the  third-order  solution  the 
average  r2  value  over  all  the  samples  was  r2  =  0.941.  For  the  first-order  solution,  the  theoretical 
solution  fit  the  measured  data  nearly  perfectly  with  an  average  r2  value  of  0.998. 
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V.  Conclusions 


The  agreement  between  the  measured  reaction  rate  data  and  the  theoretical  solution  by 
the  fitting  technique  provides  further  evidence  that  the  reaction  of  Br2  and  NO  to  form  BrNO  is 
effectively  described  by  third-order  kinetics.  Fitting  both  Ki,  which  is  directly  related  to  kf,  and 
K2,  which  is  directly  related  to  kr,  to  the  third-order  solution  produced  the  most  reliable  values  for 
kf  and  kr.  As  shown  in  Table  6  of  Section  IV,  the  value  of  kf  found  in  this  effort  is  approximately 
equal  to  previously  determined  values.  In  this  effort  kf  =  1.56  ±  0.20  x  10"38  cm6/molecule2-s  and 
in  previous  efforts  kf  ranged  from  1.18  to  1.79  x  10"38  cm6/molecule2-s,  when  including  the  error 
bounds.  The  values  of  kr  and  Keq  determined  in  this  effort  are  nearly  equal  to  the  values  found  by 
Godfrey,  but  are  not  equivalent  to  other  previously  determined  values.  In  this  effort  kr  =  2.29  ± 
0.33  x  10"21  cmVmolecule-s,  which  is  very  close  to  Godfrey’s  value  of  2.09  ±  0.55  x  10‘21 
cm3/molecule-s,  but  differs  from  earlier  values  of  2.67  and  3.71  x  10'21  cm3/molecule-s.  The 
method  used  to  solve  the  third-order  kinetics  problem  in  this  effort  is  similar  to  Godfrey’s 
technique,  which  may  partially  explain  the  close  proximity  of  the  values.  Additionally,  only 
Godfrey’s  and  this  research  accounted  for  the  loss  of  Br2  to  the  cell  walls.  The  reaction  rate 
constant  was  determined  more  accurately  by  accounting  for  this  loss. 

The  reaction  rates  determined  in  this  experiment  can  be  used  in  the  design  of  the  tunable 
laser  which  utilizes  the  NO  transition  at  5.4  microns.  The  reaction  rates  tell  how  much  BrNO  is 
present  at  equilibrium  and  can  be  used  to  determine  how  much  the  BrNO  absoiption  band  at  5.6 
microns  (v=l)  will  interfere  with  the  desired  NO  transition. 


29 


Appendix  A.  Detailed  Calculation  of  the  Kinetics  Equations 


Third-order  solution 

The  reactions  of  interest  are: 


BrNO  +  h  v  Br  +  NO 

(7) 

Br  +  BrNO  — — — >  Br2  +  NO 

(8) 

Br2  +2  NO  <  kf  >  2  BrNO 

(9) 

The  rate  of  change  of  the  concentration  of  BrNO  over  time  is  described  by 
d[BrNO\ 


dt 


■  =  -kp[BrNO\,Ip  -k\Br\[BrNO\  +2 k  f[NO\2[Br2}t -2 k,  [BrNO\2  (10) 


When  the  laser  is  shut  off  (Ip  =  0),  the  first  two  terms  go  away  and  equation  (10)  becomes 


d^^  =  2kI[NOY[Br,\-l k\BrNO\ 

at 


By  defining  the  normalized  concentration  of  BrNO  as 

..  \BrNO\ 

[ B<NO\ , 

where  the  max  value  of  x  =  1,  and  solving  for  [BrNO]t,  gives 

[BrNO],  =  x[BrNO] 


(11) 


(12) 


(19) 


Additionally, 


[NO\  =  [NO\t  +  [BrNO],,  -  [j BrNO } 
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^[NO\=[NO]eq+[BrNO]eq{\-x) 

2  [Br2  \  =  2  [Br2  ]eq  +  [l BrNO]eq  -  [ BrNO ] 

=>  \Br2  l  =  iBr2  L,  +  1  0  -  x) 


(20) 


(21) 


Substituting,  equations  (19),  (20)  and  (21)  into  equation  (11)  the  rate  of  change  of  the 
concentration  of  BrNO  over  time  becomes 


d[BrNO\ 

dt 


=  2  kf§_NO\eq  +[BrNO\eq{\-x))2{[Br2]eq  +\[BrNO\eq{\- x)\-2kr[BrNOlqx2  (32) 

v  ^ 


Using  the  chain  rule,  the  change  in  x  over  time  is 


dx  _  d\BrNO\  *  1 

dt  dt  [BrNO\e(j 


or 


dx 

dt 


2  kf[BrNOlq 


'  [NOjq 

[BrNO]eq 


J 


'  foL, 

\BrNO\q 


,  ) 

-(1-x)  -2 K[BrNO\ 


(22) 


By  letting 


[jwL 

[BrNOl, 

(23) 

K,=2  k,[BrNOl, 

(24) 

K,=2  k\BrNO\,q 

(25) 

equation  (22)  becomes 
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^  =  Ki(a  +  (l-x f 


[Br2\q  -  -2 


VN,  2 


~K2x 


(34) 


[Br2]eq/[B rN O ] eq  can  be  found  in  terms  of  a,  Ki,  and  K2  using  the  equilibrium  constant. 


7 


[flWVO] 


K"=\  [NOl,[Br,\ 


(17) 


[^L  _  kr  [BrNO]eq  _ 

2 


[BrNO]eq  kf  [NO];q  a 


(35) 


K„  = 


K. 


kf[BrNO\eq  Kj 


(36) 


Inserting  equations  (35)  and  (36)  into  equation  (34)  gives 


~T  ~  Kj(<2  +  (l-x))" 
at 


f  k2  1/,  0 


vK,cT  2 


-K2x 


J 


(26) 


The  solution  for  the  third-order  problem  is  found  by  separating  the  variables  and  integrating  both 
sides. 


\dt  =  \~ 


dx 


Kj(a  +  (l-x))" 


f  K2  1  /,  0 

vK,a  2  j 


-K2x 


(37) 


From  Mathematica,  the  solution  is 


t  _  ajfi  +/3) 
/4 


(27) 


where 


/,  =A/-K2(l  +  a)2(K2(a-l)2-2K1a3; 
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f2  =  -2(k2  -K2a2  +  Kta3 )ArcTan(K'a2{x-1+a^a2^) 
f3  =  /j (-  ln((x  - 1)2 )+  ln(Kj<22 (l  - x  +  a)2  +  2K2 (l  +  a)(l  +  x(a  - 1)+  a))) 

/4=/1(K1a3+4K2(l  +  a)) 

The  number  of  fit  parameters  in  the  solution  can  be  reduced  by  solving  Ki  in  terms  of  K2. 


K,=  2k, [BrNOl, 

(24) 

K,  =  2k,  [BrNOl, 

(25) 

[BrNO]„  =  ^ 

(38) 

^*2 kf[BrNO\„  =  K,  "K^BrNO],, 

(39) 

Using  the  relationship  in  equation  (40),  the  third-order  solution  becomes 

)/;+/;) 


a\ 


t  = 


/' 


(40) 


where 


f[  =  J-  K ,  (1  +  ■ af  (K  i  («  -  T  -  4^  2  *  Mof  ja ' ) 

f2  =  -2<K, -Km2  +(k2 

h(b  -  O2 )+  ln((K  2  *  K„  [BrNO\„  )a 2  (l  -x  +  a)2  +  2K2(l  +  «)(l  +  x(a-l)  +  a))) 

/;  =  /;  «K  2  *  [BrNOl,  )a 1  +  4K ,  (l  +  a)) 
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First-order  solution 


When  there  is  an  excess  of  NO,  the  Br2  +  2NO  2BrNO  follows  pseudo  first-order 
kinetics  and  equation  (26)  is  simplified  because  a  becomes  very  large  and  Kiot2  becomes  much 
larger  than  K2  under  these  conditions. 

f  ~T  ,  A 


^  =  K,(or  +  (l -x)y 
at 


K  2  1  /,  X 

vKj«  2  j 


-K2x 


(26) 


a  — >  oo 


K ,a2  »  K- 


K,«2  =  =  2i,[iV0]; 


[BrNO] 

Using  an  intermediate  step,  the  change  in  x  over  time  becomes 


dx  2 

—  =  K  ,a 
dt 


rw 

-K2x2  =Ky- 

f-(i 

1  -x) 

-  Kv2] 

U  J 

A2 

2 

K  xa- 

(41) 


dx  2 

—  =  K  ,a 
dt 


2  j 


Solving  the  differential  equation  gives 


(30) 


x  =  1  +  Cxe 


-t-At 


(31) 


where  A=Kj<2  and  Ci  =  -1  at  x(0)=0. 
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Appendix  B.  Kinetics  Data  Obtained  from  FTIR  spectra 


The  following  figures  show  the  kinetics  data  collected  from  the  FTIR  spectra  for  each 
sample.  The  normalized  concentration,  x,  was  plotted  as  a  function  of  time,  t.  The  normalized 
concentration  of  BrNO  was  found  by  dividing  the  area  under  the  BrNO  v=l  line  at  time=t  by  the 
area  under  the  BrNO  v=l  line  at  equilibrium.  For  each  sample  the  laser  was  turned  off  at  time=0 
sec  and  the  sample  was  allowed  to  reach  equilibrium. 


Figure  13.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  1.  Po(Br2)  =  9.07  torr,  P0(NO)  =  0.93  ton- 
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Figure  14.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  2.  P0(Br2)  =  9.08  torr,  P0(NO)  =  0.92  torr 


Figure  15.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  3.  Po(Br2)  =  9.08  torr,  P0(NO)  =  0.92  ton- 
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Figure  16.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  4.  Po(Br2)  =  9.08  torr,  P0(NO)  =  0.92  torr 


Figure  17.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  5.  Po(Br2)  =  0.25  torr,  P0(NO)  =  9.75  ton- 
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Figure  18.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  6.  Po(Br2)  =  0.30  torr,  P0(NO)  =  9.70  torr 


Figure  19.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  7.  Po(Br2)  =  0.25  torr,  P0(NO)  =  9.75  ton- 
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Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  8.  Po(Br2)  =  0.23  torr,  P0(NO)  =  9.77  torr 


Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  9.  Po(Br2)  =  0.28  torr,  P0(NO)  =  9.72  ton- 
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Figure  22.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  10.  Po(Br2)  =  0.30  torr,  P0(NO)  =  9.70  torr 


Figure  23.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  11.  Po(Br2)  =  0.30  torr,  P0(NO)  =  9.70  ton- 
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Figure  24.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  12.  Po(Br2)  =  0.25  torr,  P0(NO)  =  9.75  torr 
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Figure  25.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  13.  Po(Br2)  =  0.35  torr,  P0(NO)  =  9.65  ton- 
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Figure  26.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  14.  Po(Br2)  =  0.35  torr,  P0(NO)  =  9.65  torr 
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Figure  27.  Normalized  BrNO  concentration,  x,  plotted  as  a  function  of  time,  t,  for 
sample  15.  Po(Br2)  =  0.35  torr,  P0(NO)  =  9.65  ton- 
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Appendix  C.  Detailed  Calculations  of  kf,  kr,  Keq  values 


The  following  tables  show  the  values  that  were  used  to  determine  kf,  kr,  and  Keq  using 
first  the  third-order  solution,  equation  (27),  and  then  the  first-order  solution,  equation  (31).  First, 
considering  the  situation  where  K2  is  the  only  fit  parameter,  using  the  values  of  Po(Br2)  and 
P0(NO)  that  were  calculated  in  Section  IV  and  listed  in  Table  3,  K2  can  be  fit  to  equation  (27) 
using  Tablecurve.  K,  was  calculated  using  K2,  and  a  fixed  value  of  a  was  determined  from  [NO]o 
and  [Br2]0.  Table  9  lists  the  parameters  that  were  used  to  determine  kf  and  kr  with  only  K2  as  a  fit 
parameter.  Next,  looking  at  the  situation  where  both  Ki  and  K2  were  fit  to  the  data  using  the 
third-order  solution  and  a  was  again  fixed,  the  values  listed  in  Table  10  were  used  to  determine 
kf,  kr  and  Keq.  Finally,  looking  at  the  situation  where  A=Kia2  was  fit  to  the  collected  kinetics  data 
using  equation  (31),  the  values  listed  in  Table  1 1  were  used  to  calculate  kf  and  kr. 


Table  9.  Values  used  to  determine  kfand  kr  for  case  where  K2  was  the  only  fit  parameter 


f(NO)0 

(unitless) 

f(Br,)0 

(unitless) 

neq 

(torr) 

a 

(unitless) 

[BrNO]Cq 

(molecule/cm3) 

[NO],.q 

(molecule/cm3) 

KeqBrNOc, 

(unitless) 

K, 

(sec"1) 

k2 

(sec'1) 

K2  Std 
Error 
(sec'1) 

r 

Px  =  10.0  torr,  T=293  K,  Keq  =  6.7 IE- 18  cm3/molecule 

0.0925 

0.9075 

0.2912 

0.717 

1.78E+16 

1.27E+16 

0.119 

8.79E-6 

7.38E-5 

6.80E-7 

0.985 

0.0920 

0.9080 

0.2913 

0.717 

1.77E+16 

1.27E+16 

0.118 

8.72E-6 

7.36E-5 

5.66E-7 

0.983 

0.0922 

0.9078 

0.2913 

0.717 

1.77E+16 

1.27E+16 

0.119 

8.40E-6 

7.08E-5 

3.27E-7 

0.996 

0.0920 

0.9080 

0.2913 

0.717 

1.77E+16 

1.27E+16 

0.118 

7.71E-6 

6.51E-5 

4.52E-7 

0.989 

0.9750 

0.0250 

0.2381 

19.477 

1.57E+16 

3.06E+17 

0.105 

6.95E-6 

6.60E-5 

1.53E-6 

0.928 

0.9700 

0.0300 

0.2827 

16.158 

1.86E+16 

3.01E+17 

0.125 

1.02E-5 

8.17E-5 

2.79E-6 

0.863 

0.9750 

0.0250 

0.2381 

19.477 

1.57E+16 

3.06E+17 

0.105 

7.90E-6 

7.51E-5 

1.81E-6 

0.920 

0.9770 

0.0230 

0.2199 

21.216 

1.45E+16 

3.07E+17 

0.097 

6.79E-6 

6.99E-5 

1.88E-6 

0.908 

0.9720 

0.0280 

0.2650 

17.343 

1.75E+16 

3.03E+17 

0.117 

8.48E-6 

7.24E-5 

1.39E-6 

0.952 

0.9700 

0.0300 

0.2827 

16.158 

1.86E+16 

3.01E+17 

0.125 

1.11E-5 

8.86E-5 

1.22E-6 

0.959 

0.9700 

0.0300 

0.2827 

16.158 

1.86E+16 

3.01E+17 

0.125 

1.11E-5 

8.91E-5 

9.79E-7 

0.949 

0.9750 

0.0250 

0.2381 

19.477 

1.57E+16 

3.06E+17 

0.105 

8.54E-6 

8.11E-5 

1.34E-6 

0.942 

0.9650 

0.0350 

0.3262 

13.791 

2.15E+16 

2.97E+17 

0.144 

1.56E-5 

1.08E-4 

2.54E-6 

0.908 

0.9650 

0.0350 

0.3262 

13.791 

2.15E+16 

2.97E+17 

0.144 

1.74E-5 

1.21E-4 

2.63E-6 

0.930 

0.9650 

0.0350 

0.3262 

13.791 

2.15E+16 

2.97E+17 

0.144 

1.70E-5 

1.18E-4 

2.98E-6 

0.893 
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Table  10.  Values  used  to  determine  kf  kr  and  Keq  for  case  where  and  K2  were  the  fit 
parameters 


f(NO)o 

(unitless) 

f(Br2)0 

(unitless) 

neq 

(torr) 

a 

(unitless) 

[BrNO]cq 

(molecule/ 

cm3) 

[NO]cq 

(molecule/ 

cm3) 

K, 

(sec'1) 

K,  Std 
Error 
(sec'1) 

k2 

(sec'1) 

K2  Std 
Error 
(sec*1) 

r 

PT  =  10.0  torr,  T=293  K,  Kcq  =  6.7 IE- 18  cm3/molecule 

0.0925 

0.9075 

0.2912 

0.717 

1.78E+16 

1.27E+16 

9.34E-6 

4.93E-8 

7.87E-5 

4.23E-7 

0.997 

0.0920 

0.9080 

0.2913 

0.717 

1.77E+16 

1.27E+16 

8.56E-6 

9.36E-8 

7.22E-5 

8.07E-7 

0.984 

0.0922 

0.9078 

0.2913 

0.717 

1.77E+16 

1.27E+16 

8.30E-6 

4.77E-8 

6.99E-5 

4.08E-7 

0.997 

0.0920 

0.9080 

0.2913 

0.717 

1.77E+16 

1.27E+16 

7.38E-6 

5.64E-8 

6.23E-5 

4.83E-7 

0.993 

0.9750 

0.0250 

0.2381 

19.477 

1.57E+16 

3.06E+17 

6.87E-6 

2.74E-7 

6.75E-5 

4.71E-6 

0.928 

0.9700 

0.0300 

0.2827 

16.158 

1.86E+16 

3.01E+17 

1.07E-5 

6.16E-7 

7.49E-5 

6.65E-6 

0.865 

0.9750 

0.0250 

0.2381 

19.477 

1.57E+16 

3.06E+17 

7.95E-6 

3.41E-7 

7.41E-5 

6.00E-6 

0.920 

0.9770 

0.0230 

0.2199 

21.216 

1.45E+16 

3.07E+17 

6.84E-6 

3.15E-7 

6.88E-5 

5.86E-6 

0.908 

0.9720 

0.0280 

0.2650 

17.343 

1.75E+16 

3.03E+17 

8.59E-6 

2.76E-7 

7.06E-5 

3.57E-6 

0.952 

0.9700 

0.0300 

0.2827 

16.158 

1.86E+16 

3.01E+17 

1.13E-5 

2.28E-7 

8.56E-5 

2.28E-6 

0.960 

0.9700 

0.0300 

0.2827 

16.158 

1.86E+16 

3.01E+17 

1.20E-5 

3.23E-7 

8.01E-5 

2.94E-6 

0.941 

0.9750 

0.0250 

0.2381 

19.477 

1.57E+16 

3.06E+17 

8.11E-6 

2.27E-7 

8.87E-5 

3.99E-6 

0.944 

0.9650 

0.0350 

0.3262 

13.791 

2.15E+16 

2.97E+17 

1.54E-5 

6.42E-7 

1.13E-4 

9.05E-6 

0.908 

0.9650 

0.0350 

0.3262 

13.791 

2.15E+16 

2.97E+17 

1.73E-5 

6.25E-7 

1.23E-4 

8.34E-6 

0.930 

0.9650 

0.0350 

0.3262 

13.791 

2.15E+16 

2.97E+17 

1.72E-5 

7.70E-7 

1.14E-4 

1.04E-5 

0.893 

Table  11.  Values  used  to  determine  kfand  kr  for  case  where  K,a2  was  the  fit  parameter 


^NOo 

(unitless) 

fl3r2o 

(unitless) 

neq 

(torr) 

[BrNO]cl| 

(molecule/ 

cm3) 

[NO]eq 

(molecule/ 

cm3) 

K,a2 

(sec'1) 

K,a2  Std 
Error 
(sec'1) 

a 

(unitless) 

K1 

(sec'1) 

K2 

(sec'1) 

r2 

PT  =  10.0  torr,  T=293  K,  Keq  =  6.7 IE- 18  cm3/molecule 

0.9750 

0.0250 

0.2381 

1.57E+16 

3.06E+17 

2.66E-3 

3.54E-5 

19.477 

7.00E-6 

6.65E-5 

0.997 

0.9700 

0.0300 

0.2827 

1.86E+16 

3.01E+17 

2.94E-3 

3.51E-5 

16.158 

1.13E-5 

9.02E-5 

0.998 

0.9750 

0.0250 

0.2381 

1.57E+16 

3.06E+17 

2.65E-3 

3.21E-5 

19.477 

6.97E-6 

6.63E-5 

0.997 

0.9770 

0.0230 

0.2199 

1.45E+16 

3.07E+17 

2.46E-3 

3.19E-5 

21.216 

5.47E-6 

5.63E-5 

0.997 

0.9720 

0.0280 

0.2650 

1.75E+16 

3.03E+17 

2.67E-3 

2.50E-5 

17.343 

8.86E-6 

7.57E-5 

0.998 

0.9700 

0.0300 

0.2827 

1.86E+16 

3.01E+17 

3.03E-3 

2.00E-5 

16.158 

1.16E-5 

9.28E-5 

0.998 

0.9700 

0.0300 

0.2827 

1.86E+16 

3.01E+17 

2.96E-3 

1.85E-5 

16.158 

1.13E-5 

9.06E-5 

0.999 

0.9750 

0.0250 

0.2381 

1.57E+16 

3.06E+17 

2.99E-3 

2.00E-5 

19.477 

7.88E-6 

7.49E-5 

0.999 

0.9650 

0.0350 

0.3262 

2.15E+16 

2.97E+17 

5.28E-3 

7.00E-5 

13.791 

2.77E-5 

1.92E-4 

0.997 

0.9650 

0.0350 

0.3262 

2.15E+16 

2.97E+17 

5.14E-3 

7.40E-5 

13.791 

2.70E-5 

1.87E-4 

0.997 

0.9650 

0.0350 

0.3262 

2.15E+16 

2.97E+17 

5.21E-3 

6.80E-5 

13.791 

2.74E-5 

1.90E-4 

0.997 
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